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Abstract 

The quantum Josephson Hamiltonian of two weakly linked Bose-Einstein 

condensates is written in an overcomplete phase representation, thus avoiding 

the problem of defining a Hermitian phase operator. We discuss the limit of 

validity of the standard, non-rigorous Mathieu equation, due to the onset 

of a higher order cos 20 term in the Josephson potential, and also to the 

overcompleteness of the representation (the phase (p being the relative phase 

between the two condensates). We thereby unify the Boson Hubbard and 

Quantum Phase models. 
PACS: 74.50.+r,03.75.Fi 
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Recent developments on the engineering of Bose-Einstein condensates (EEC's) |T|,0] are 
suggesting new scenarios for searching macroscopic quantum coherence phenomena. Most of 
the dynamical regimes investigated experimentally have been understood within the Gross- 
Pitaevskii (GP) framework 0. GP is essentially a "classical" theory, since it approximates 
the quantum bosonic field as a complex order parameter the new challenge is the 
systematic study of regimes where quantum fluctuations cannot be ignored. 

Quantum effects play an important role in a "mesoscopic" Josephson junction P]. The 
Josephson effects (JE) are a paradigm of phase coherence in superfluid/superconductive 
systems In the BEG context, JE have been associated with coherent collective oscilla- 
tions between two weakly-linked condensates, the weak link being created by the potential 
barrier in a double- well trap [0-0, or by a quasi- resonant external elettromagnetic field, 
which induces oscillations between two trapping hyperfine levels Although bosons 

are neutral, and external circuits are obviously absent, it is also possible to study the close 
analougs of the the "ac" and "dc" effects observed in superconducting Josephson junctions 
|16| . Quite recently, an array of "mesoscopic" weakly linked condensates has been created 
in [|^, each trap, located at the antinodes of an optical standing wave, containing ~ 1000 
condensate atoms. 

The classical Josephson Hamiltonian (GJH) has been casted in term of a pendulum-like 
equation ||7|,p!7|: 



i^d = yn^-^^iVVl^cosc/. (1) 

with the relative number of condensate atoms between the two bulk condensate n = ^{ni — 
playing the role of the momentum, and the relative phase = 0i — 02 being the angle 



respect to the horizontal axis |18 



The "charging energy" and the "Josephson coupling energy" Ej can be calculated as 
overlap integrals |7|,p!6|: 



Ec = 2g J dr $i(r)^ = 2^/ ^^(r)^ (2a) 
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Ej = -N J dr $i(r) 



$2(r) 



(2b) 



with the one-body wave functions $i(r),$2('") locahzed in the trap 1,2, respectively, and 
/ dr $i(r)*$2(r) = 0, / rfr 2(^^)1 = 1, fi' = being the scattering length and m the 

atomic mass; is the total number of atoms |T^. 

The quantum Josephson Hamiltonian (QJH), on the other hand, has been studied within 



two different models which seem in conflict with each other |£0| . This conflict originates, at 
its heart, from a fundamental, yet unsolved, problem: the existence and the meaning of a 



quantum phase operator |21 



In the "boson mode" representation, QJH is written in terms of creation and destruction 
boson operators [P,p!^, p!^ , p!5| ] : 



E E 
H = — (cij'^aj'^aidi + 0,^0,^0,202) — j;^{d^CL2 + o^Oi) 



(3) 



where 0^2, (^'1,2) creates (destroys) a particle in the trap 1, 2, respectively. 

In the "phase" representation, on the other hand, the relevant quantum observables are 
the difference of phases and number of atoms between the two condensates in each trap 
pn] , |22|j23| , the Hamiltonian being written in term of a "quantum pendulum" (Mathieu) 
equation: 

E., d' 



H 



Ej cos ( 



(4) 



with fi^ = (dj^di — 0^0,2)^ ~ ^ which cannot, as we will discuss shortly, be 

unambiguously defined in terms of creation/destruction operators. Eq. (^ acts on a 27r- 
periodic "wave-function" V(0) = V'(0 + 2vr), with ^ jl-^dcp ip^cj))* 4'{4>) = 1 and expectation 
values {A) = ^ jl-^dcj) A{(j),t) tp{(f)). This framework, generalized to include damping 

effects and external circuits, is rather standard in the SJJ literature ||2^ . 

Generally speaking, the Hamiltonian Eq. (|^) can be derived microscopically from the 
quantum field equation: 



H 



j dr if' 



— V" + Kxt + - 
2m 2 



(5) 



writing the boson field in the "two-mode approximation" as: ^ ~ o-i + $2('") 0,2- 

The quantum pendulum Eq. on the other hand, can be retrieved quantizing the 
classical Josephson Hamiltonian Eq. (|l|) by replacing the classically conjugate phase/atom 
number with the corresponding, non-commuting, operators . 



On this grounds we can consider Eq. as a phenomenological, while Eq. has a 
clearer "microscopical" foundation: it is generally believed that, for some limit, these two 
Josephson Hamiltonians describe essentially the same physics, yet, their exact relation has 
not been exploited in the literature. A way to analyze, a posteriori, the relation between 
the two Hamiltonians Eq. (^) and Eq. (^, is to identify: 

a, = e*ril/' (6) 

with 0j and = afai Hermitian phase and number operators acting in the trap i. The 
assumed Hermiticy of 0j implies the unitarity of exp{i(j)i). Number and phase, thus, are 
treated as conjugate observables: 

[kni]=t (?) (7) 

which implies, in the phase-representation, rij = ~i-^,(f)i = (pi. This approach is known 
to be incorrect. The commutator gives rise to inconsistencies when its matrix elements are 
calculated in a number-state basis. Even worse, the exponential operator exp{i(f)) is not 



unitary and so does not define an Hermitian ||2 1|| . 

Despite such problems, the Josephson Hamiltonian, Eq. (^), is considered as the starting 
point of most analysis in SJJ and BJJ, the implicit caveat being that the commutator 
relation Eq. (0) is approximately correct for systems with a large number of "condensate" 
Cooper-pairs/atoms. More precisely, the relation between Eq. (|^) and Eq. is considered 
only in a semiclassical context, replacing di^2 ^ ^^■1,2^^'^^'^ , ^1,2 being c-numbers, in the term 
proportional to the Josephson coupling energy Ej in Eq. (^). 

In this letter we develops a consistent procedure to rewrite the quantum Josephson 
Hamiltonian Eq. @ in the phase representation, overcoming the problem of defining an 



Hermitian phase operator. The resulting equation differs in several ways from the Mathieu 
equation Eq. which is recovered in a specific limit. Although we consider, for simplicity, 
the case of a two-sites Josephson junction, our analysis clarifies the relation between the 
Boson Hubbard p^j25|] (which can be seen as the n-sites generalization of Eq. (|^)) and the 



Quantum Phase models (n-sites generalization of Eq. (Q)) governing the dynamics of a 
EEC's array as that created in [0. 

An arbitrary state spanning the two-dimensional Hilbert space, in which is defined the 
QJH Eq. (H), can be expanded as: 

1 r 

rf0lrf02 / (01, 02)101, 02) (8) 



with the (un-normalized) Bargmann state 

OO -| 

101,02)= E e^('<^^+-^^)^=|0|m) (9) 

with |m), |/) atom number eigenstates of trap 1, 2 respectively. 

In the Bargmann space, the action of the boson operators on the state vector Eq. (^ 
reduces to differential operators acting on /(0i,02): 

/(0,) a, e*^' /(0,) (10a) 

/(0.)a+^z^[e-^7(0.)] (10b) 
d<pi 

as can be seen after an integration by parts. In particolar, the atom number operator 
assumes the familiar form: /(0j) fii = f{(f)i) afai ""^^ /(0j)- 
We consider as fixed the total number of atoms, and we can write: 

/(0,0+) = e-^^^+^(0) (11) 

with 0+ = |(0i + 02), = 01 — 02, iV = ni + n2, n = ^{rii — 712), and n being the total 
and the relative number of atoms. ■0(0) is an arbitrary, 27r-periodic function normalized to 
unity: ^ /Z!^ (i0|'0(0)p = 1. We can integrate over 0+ the Eq. (||): 

1 



2tt 



rf0^(0)|0) (12) 



with 

|^)= f; -j=^2M==\n) = m(ate'Uate-''rf\0) (13) 
that can be considered, in the large N hmit, as a pure phase state [^. The state Eq. (13) 



is proportional (with = 0) to the exact ground state of Eq. (0) in the non-interacting limit 
{Ec = 0). Here we consider the states |0) as an overcomplete base to expand the solution 
of the full Josephson Hamiltonian Eq. (^. The scalar product of two phase states is given 
by: 

i^l^) = 7^ r dcj) ^Piey V^(0) (^^10) = 1 (14a) 



(27r)2 J- 

In the limit of large number of atoms, {0\(j)) becomes proportional (in the interval —it < 
(f), 6 < it) to a delta function: 

lim (5^(0 -6) = 6{(j) - 6) (15) 
The action of the Josephson Hamiltonian Eq. (|^) on the state vector Eq. (0) gives: 

H\i,) = {H, + Hjm (16a) 

m) = jy^ \^) [-2^ + - N)] m (16b) 

Hj\i^) = 10) [(iV + 2)cos0 + 2sin0 ^] V^(0) (16c) 

Projecting out in the phase space: 

(^1^1^) = i- £^ d<f> F ^(0) (5|0) (17a) 
^V'(0) = [-Y^-i^./(l + |)cos0 -2^ sin0^]^(0) (17b) 

where a constant energy shift has been omitted. The spectrum of H is given by the eigenvalue 



equation for this F, which is known as the Ince equation p8|. Notice that, because of the 
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overcompleteness of the phase-states Eq. (|T3|), Hj in Eq. (|T6| ,c) is Hermitian, even though 
the differential operator sin0 ^ is not. 

Thus the time-dependent Schrodinger equation in our phase representation, 

i^^m) - mH,+ = Q , (is) 

can be put in dimensionless form (by defining 7 = and rescahng t t^) as 

£j<j)6N{<f)-e)[t- + — + ^{N + 2)cos<j) +27 sin0— ] V'(0) = O (19) 

We may obtain a more transparent form of (p!9D by defining ip{4>) = e'^^°^'^'^{(j)). This 
ehminates the first derivative and, dropping a constant, maps F onto a manifestly Hermitian 
operator acting on ^f. The result is that to satisfy Eq. (|^) it is sufficient to satisfy 

= [-^ - ^{N + 1) cos - -7' cos 20]vl> . (20) 

Setting = gives the three-term Hill equation |2 



Eq. (|20| ) is of the same usefully simple form as the QPM. One can apply to it all the 
intuition and experience, and all the calculational tools, that are relevant to the one-particle 
Schrodinger equation with a periodic potential. Solving it yields the energy spectrum of 
the BHM, and thus any dynamical time scales (such as for decay of the possible metastable 
states by tunneling). But since the equation is based on an overcomplete representation, 
expectation values and matrix elements must be computed using the nonstandard inner 
product with the kernel 6n, and with the extra e'^^"^'^ factors; hence the interpretation of 
the eigenstates can be somewhat subtle. Since 6j\f, with >> 1, behaves as a delta function 
for test functions varying slowly on the scale N~^^'^, for the lower energy states the inner 
product is essentially the standard one, and Eq. (|20| ) can be used without any extra thought. 
For the higher energy states, however, the interpretations of "^{(p) and Eq. (^) depend on 
7- 

We can distinguish three regimes. For 7 > 0{N), the cos 20 term in (^) is a significant 
correction to the purely cos potential of the QPM; and for 7 > N/2 it makes a second local 



minimum at = vr. This agrees with the GPE result that for 7 > N/2 there are metastable 
TT- oscillations of zero average n. But the trapping of probability near = vr is not trivial, 
because e^^"^'^ factor tends strongly to suppress it. This suppression is countered for high 
energy states, however, by the fact that 5n effectively eliminates all Fourier components e^^"^ 
with \k\> N/2 in ip{(t)). (This is as it should be, since from our derivation it is clear that such 
components are unphysical.) There are thus high energy solutions to Eq.(pO|) whose physical 
components are localized around (p = n, because in the deeper well around = their WKB 
frequency has \k\ > N/2. In the intermediate range of energies, however, the implications 
of 6n are more complicated. Eqn. (^) may still be treated with semiclassical methods, as 
will be discussed elsewhere. It is qualitatively clear that both the second minimum of the 
potential, and the failure of quantum motion to be strictly confined to it, will be important; 
and so we must consider (^) as giving corrections to both the QPM and the GPE. 

For 1 < 7 << A^, the cos 20 term is insignificant, and so the eigenf unctions of Eq. (^) are 
essentially the Mathieu functions of the QPM. The nonlocality of 6n is of little qualitative 
significance, except for the highest energy states. These are well approximated by the 
WKB method, and their local wave number is of the form k + 7"^^^ cos0. For k of order 
N/2, this implies that solutions to (|2^) of the 'running' type, with no localized phase, 
actually have their physical component localized around = tt. This supports the GPE 
in what is its only significant discrepancy with the QPM in this range of 7, for the GPE 
also predicts high energy vr-states with nonzero average n [0. (We neglect the fact that 
the actual quantum eigenstates must have zero average n, because they are even and odd 
'Schrodinger's Cat' superpositions of the states with large positive and negative n. But the 
splitting is exponentially small in A^.) 

For 7 < 1, the situation is similar, except that the modulation of the local wave number 
with is too small for the elimination of unphysical Fourier components to produce a 
localized state. It simply eliminates all WKB states with \k\ > N/2, and the result is that 
the physical spectrum of Eq. ( pOD is essentially Mathieu functions, up to a cut-off. The vr- 
states of the GPE are absent and the QPM is vindicated, except for the welcome truncation 
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of the spectrum, in accordance with the fact that the Hilbert space must indeed be only 
(A^ + l)-dimensional. 

In Fig. we draw the probabihty in phase, P{(f)), in an orthonormal representation 
of (f). This illustrate the interpretation of our nonorthonormal phase representation with 
unphysical components projected out. Considering = 199 particles and 7 = all 200 
energy eigenstates are found numerically in the number basis, and ranked in order of energy. 
The absolute square of the Fourier transform of the number amplitudes is plotted for each 
state; since all P{4>) are even, only < < vr is shown. This approach is not possible 
for large A^, but it illustrates and confirms our discussion of the physical interpretation of 
our overcomplete representation, which is useful for all A^. For this value of 7, Josephson- 
localized, running, and vr-localized states all occur. 

To conclude, a few words are in order about the Rabi limit 7 — > 00. In this case 
the potential of Eq. (20) becomes a pair of strong harmonic potentials centred at = 
and = TT, both having the frequency 'jEc = 2Ej/N. The equally spaced harmonic 
oscillator eigenstates of these potentials, up to the cut-off described earlier, reproduce the 
Rabi spectrum of the Hamiltonian (3) in the limit —* 0. 

This work has been partially supported by the Cofinanziamento MURST. 
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FIGURES 

FIG. 1. The module square of the energy eigenstates P{(f)), as a function of an orthonormal 

representation of (p, with N = 199 and 7 = 
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